Electromagnetic response of layered superconductors with broken lattice inversion 

symmetry 
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We investigate the macroscopic effects of charge density waves (CDW) and superconductivity in 
layered superconducting systems with broken lattice inversion symmetry (allowing for piezoelectric- 
ity) such as two dimensional (2D) transition metal dichalcogenides (TMD). We work with the low 
temperature time dependent Ginzburg-Landau theory and study the coupling of lattice distortions 
and low energy CDW collective modes to the superconducting order parameter in the presence of 
electromagnetic fields. We show that superconductivity and piezoelectricity can coexist in these sin- 
gular metals. Furthermore, our study indicates the nature of the quantum phase transition between 
a commensurate CDW phase and the stripe phase that has been observed as a function of applied 
pressure. 
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I. INTRODUCTION 

The quasi-two-dimensional (2D) transition metal 
dichalcogenides (TMD)i 2H-NbS 2 , 2H-NbSe 2 , 2H-TaS 2 
and 2H-TaSe 2 are layered compounds where supercon- 
ductivity coexists with a charge density wave (CDW) 
stated. It has been documented experimentally that 
these materials present anomalous properties such as a 
decrease of the resistivity and decoupling between pla- 
nar and c-axis transport in the CDW phase 3 , anomalous 
impurity effects on the superconductivity and non-linear 
Hall effect, stripe phases^, and more recently angle re- 
solved photoemission experiments (ARPES) have shown 
a quasi-particle lifetime diverging linearly with the en- 
ergy close to the Fermi surface^ in contrast to ordinary 
Fermi liquids where the lifetime diverges quadratically 
with energy. Some of these properties are similar to 
the ones observed in high temperature superconductors 
(HTc)i. Contrary to HTc, however, TMD are extremely 
clean systems and therefore the anomalous behavior is 
undoubtedly intrinsic and sample independent. There- 
fore the understanding of the anomalous metallic behav- 
ior in these 2D materials may provide clues for the under- 
standing of the anomalous physics that occurs in a broad 
class of 2D systems where superconductivity occurs. 

One of us (AHCN), has proposed recently an uni- 
fied picture of the CDW and superconducting transi- 
tions that can explain some of the anomalies observed 
experimentally^. The theory proposes that the mech- 
anism of superconductivity is related to the pairing of 
the CDW elementary excitations mediated by acoustic 
phonons via a piezoelectric coupling. These excitations 
are Dirac electrons situated in the nodes of a CDW gap, 
resembling the Fermi surface zone of graphite. In TMD, 
unlike the case of graphite, the lattice inversion sym- 
metry is broken in the CDW phase. The breaking of 
the inversion symmetry in TMD, as documented in neu- 



tron scattering experiments 9 , allows for a piezoelectric 
coupling between charge carriers and phonons 10 . Since 
CDW formation is usually related to nested Fermi sur- 
faces in ID systems, weak 2D nesting combined with 
strong variations in the electron phonon coupling due to 
the tight-binding nature of the electronic orbitals, can be 
responsible for the origin of the nodal CDW order param- 
eter. This model is able to correctly explain some of the 
anomalous properties of TMD such as the energy depen- 
dence of the quasi-particle lifetime (given by the imagi- 
nary part of the Dirac fermion self-energy), the critical 
dependence of the superconducting transition with the 
lattice parameters, the reduction of the resistivity, and 
the metallic behavior in the CDW phase. 

However, a criticism has been raised to this theory 
based on the fact that piezoelectricity normally occurs 
in insulating systems since metals can screen the inter- 
nal electric fieldsii. However, Dirac fermions in 2D have 
a vanishing density of states at the Fermi energy and 
therefore do not screen electric fields^ allowing for the 
existence of a metallic state (that is, a state with gapless 
fermionic excitations) and piezoelectricity. Nevertheless, 
the possibility of coexistence of piezoelectricity and su- 
perconductivity is indeed surprising. In such a system it 
would be possible to generate super-currents by simply 
squeezing the sample. In this work we show that there is 
no contradiction between the existence of piezoelectricity 
and superconductivity in materials where the low-lying 
excitations are Dirac fermions. We investigate the cur- 
rent and charge fluctuations that would arise in these 
systems under the view point of the collective modes. 

Using path integrals we derive a semi-classical action 
that describes the coupling between plasmons (responsi- 
ble for screening) and acoustic phonons to the supercon- 
ducting order parameter via the piezoelectric coupling. 
Various well-known regimes are described by this action: 
(1) in the case of normal electrons without supercon- 
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ductivity we recover the well-known results for collective 
modes and screening in 2D and 3Di 3 «; (2) in the case of 
Dirac fermions without superconductivity we re-obtain 
the results for collective modes and screening in semi- 
metals like graphite-i; (3) in the absence of piezoelectric 
coupling, we recover the behavior of an ordinary type II 
superconductor 15 . However, when we allow the piezo- 
electric coupling with Dirac fermions and superconduc- 
tivity, new effects appear in the electromagnetic response. 
Moreover, our results shed light on the origin of the quan- 
tum phase transition between the commensurate CDW 
phase and the stripe phase in TMD. We also investigate 
the collective modes that appear when combined with 
low lying energy bulk plasmons. 

The article is organized in the following manner: in 
section II we discuss the plasmons in layered nodal liquids 
and specifically in TMD; in section III we develop the 
semi-classical calculation in general grounds and apply it 
to layered compounds like TMD in section IV. We have 
added an appendix where the details of the calculation 
are presented. 



II. NODAL LIQUID PLASMONS 

The field theory of nodal liquids^ is built under the 
basic idea that the nodes of a CDW order parame- 
ter lead to two distinct subsystems which correspond 
to the two components of the Dirac fermion spinor 
= (ip' +a ,ip_ a ), with +,— indexing respectively the 
fermionic particles and anti-particles (holes). The non- 
interacting low energy Hamiltonian that describes the 
elementary excitations inside the Dirac cone is: 



/d 3 k 
j^^$i(k) h(v F a x k x + VA<Tyk y ) <f> a (k) , 

w 

where k is the momentum, a x , v are Pauli matrices that 
act in the particle-anti-particle subspace, vf and va are 
the anisotropic velocities perpendicular and parallel re- 
spectively to the Fermi surface, and BZ is the first Bril- 
louin zone of size 2?r /d along the k z axis (d the inter-plane 
distance). When the chemical potential, /x, intercepts the 
Dirac point, the Hamiltonian leads to the zero order 
polarization function^ at T = 0, whose complex conju- 
gate is given by: 
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where u> is the frequency, q = q x + {va,/vf)% the 
anisotropic in-plane momentum. In the absence of hop- 
ping between the planes or interactions with the lattice, 
the collective excitations of the Dirac fermions are due 
to the 3D Coulomb interaction, Vq, between carriers in 
different planes: Vo(q,k z ) = 2irde 2 /(e q)S(q,k z ), with a 
structure factor S(q, k z ) = smh(qd)/[cosh(qd) — cos(k z d)] 
identical to the layered electron gas (LEG) casein (e is 



the electron charge, and eo the dielectric constant). In 
the random phase approximation (RPA) the electric sus- 
ceptibility is given by: e(w,q, k z ) = l-V Q (q, fc 2 )n°(cj,q). 
At the points where the susceptibility vanishes, that is, 
e(w p (q, k z ), q, k z ) = 0, one obtains the plasmon disper- 
sion, w p (q, k z ). From JSJl, we notice that no plasmon 
modes are allowed for Dirac fermions. 

The displacement of the chemical potential from the 
Dirac point generates a pocket around the nodes which 
drastically changes this picture. The density of states 
becomes finite at the Fermi surface with Fermi momen- 
tum k F , and Fermi energy E* F = HkpVF, giving rise to 
intra-band excitations in the cone. In this case we re- 
cover the optic plasmon 0J p (q) — ^Juj 2 {q/q) 2 + (voq) 2 for 
the k z = bulk mode, where vq is the plasmon speed, 
and the 2D acoustic mode Lu p (q) cx y/q 2 /q for the rest of 
the plasmon band < k z < Tr/d. The intra-band con- 
tribution is satisfactorily understood in the context of 
doped graphite^ where for vf<1 < uj < VF{2k F — q) there 
is gap in the particle-hole continuum of the inter-band 
excitations (see fig.l) defined by the imaginary part of 
@. The optical plasmon is free of Landau damping in 
the long wavelength limit and its energy is of the order 
of E F . Notice that in the anisotropic case (vf ^ va) the 
gap in the plasmon spectrum depends on the direction 
around the pocket Fermi surface. 




FIG. 1: Schematic drawing of the k z = bulk plasmon in 
an isotropic pocket with Fermi momentum k* F . v — uj/vf 
and q is the in-plane transfer momentum. The shaded areas 
correspond to the particle-hole continuum due to intra-band 
and inter-band excitations. 

If we also add the piezoelectric electron-phonon 
coupling 8 



Hep = iY.J d'^W 4(x)$ CT (x) 
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due to acoustic phonons (with phonon field </>(x)) with 
energy hujq into the pocket excitations, the RPA elec- 
tric susceptibility acquires a correction: e(w,q, k z ) = 
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1 - [V (q,k z ) + ( 7 2 //l)Z)V,q)]n (u;,q,M) where 



13 



D°(cj,q) 



- (w q - ill) 2 



is the phonon propagator that affects very little the op- 
tical bulk plasmon in the q — ■> limit. 

Besides doping, coherent inter-layer hopping is known 
to drive a dimensional crossover in direction to a 3D sys- 
tem, leading as well to a pocket formation in a nodal 
liquid—. The energy of the pocket is of the same order 
of the inter-layer hopping energy, which in TMD NbSe2 
and NbS2 (we will drop the 2H prefix from now on) were 
calculated in ~ 0.1 eV 2 ^. We notice that inter-layer in- 
teractions were not taken into account into those band 
calculations and that this energy could be considerably 
smaller. More experimental studies are required to in- 
vestigate the nature of the low energy collective modes 
in these materials. However, we can consider the plas- 
mon modes in real materials to have a finite energy, like 
in ordinary metals. Nevertheless, the plasmon frequency 
will be smaller than in ordinary metals because of the 
low density of states in the system. 



III. SEMI-CLASSICAL DYNAMICS 

We consider the problem of Dirac fermions described 
by flj coupled to lattice vibrations described by the 
Hamiltonian 
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where n labels the lattice site B? n , P* is the lattice mo- 
mentum operator with i = x,y, z (these operators are 
canonically conjugated, namely, [P^,R J rn ] = ifr5ij5 nrn ), 
K nm is the coupling matrix. In what follows we will 
also consider the coupling of the Dirac fermions and lat- 
tice degrees of freedom to a classical electromagnetic field 
described by a vector potential field, A(r), and electro- 
magnetic field tensor, F» v = d^A v - d v A» (d» = d/dx^ 
with fi, v — 0, i) with electromagnetic energy density: 
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where Foj — ~Ej is the electric field strength (we use 
the metric g^ >u — (1, — 1, — 1, — l)) 21 . In a piezoelectric 
the electric field couples to the lattice distortion, X l n = 
R l n — R l n+ i, via the piezoelectric tensor, A*-, by: 



(6) 



where repeated indices are to be summed. The coupling 
between Dirac fermions and lattice vibrations is given by 
@. The full problem can be written in path integral 
form for the generating functional of the problem: 



dt / dr£[^,X,A M ] \ ,(7) 



Z = I DipDip I L>Xexp<{ - 



where C is the Lagrangian (real time) of the problem (-0, 
tp are Grassmann variables). 

Let us consider first the case of the normal (non- 
superconducting) state. We are only interested in the 
long wavelength physics of the electronic problem which, 
as argued in the previous section, is described by the plas- 
mon mode. The plasmon mode can be separated from 
the particle-hole continuum using the formalism devel- 
oped by Bohm and Pines 2 ^ and the Lagrangian of the 
problem reduces to: 
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where Q w Air/k 2 , with k c the 3D electron screening wave 
vector (for an isotropic metal k c — y^Snne 2 / (2Ep) where 
n is the 3D electron density and Ep the Fermi energy 2 ^), 
k is the lattice mass density, v p h is the sound velocity, 
j$/c = pn is the density of normal electrons and 
is the normal Ohmic current density (in our notation, 
d° = c~ 1 d t is the time derivative normalized by the speed 
of light c) . 

The collective properties of the electrons in the nor- 
mal phase are the same of the electrons in the supercon- 
ducting phase, since the opening of the superconducting 
gap does not affect the real part of the electronic polar- 
ization function, n(q), in the long wavelength limit 23 . 
In fact, within RPA one can easily show in the con- 
text of BCS theory that Ren(q) remains essentially un- 
changed until the superconducting gap is of the order of 
the Fermi energy, when the BCS pairing approximation 
becomes invalid 2 *. In other words, for small q the plas- 
mon is not sensitive to the superconductor phase transi- 
tion. Thus, normal and superconducting electrons screen 
electric fields exactly in the same way and therefore enter 
at equal footing in the Lagrangian. Thus, irrespective of 
the phase, normal or superconducting, p^ in JSJ can be 
replaced by the total electron density p = p^ + p s , where 
p s is the density of superconducting electrons (p s = 
in the normal state). Normal currents, however, are due 
to the quasi-particle excitations while the supercurrents, 
j/, are due to the ground state of the condensate 2 ^. 
Thus, in the superconducting state the vector potential 
only couples to the normal currents. At finite tempera- 
ture, T, the number of normal and superconducting elec- 
trons is not conserved. For simplicity, we have chosen to 
work at T = 0, where all the electrons are in the conden- 
sate and the total current is — jf". 

In the superconducting phase we can introduce the 
superconducting order parameter, 'J, via a standard 
Hubbard-Stratanovich transformation and trace over the 
electrons^. The generating functional reads: 



Z = I D^D^ I £)Xexp<{ - / dt / dr(£ N [p,X, A] 



£ GL [*,A])} 



(9) 
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where Cgl is the Ginzburg-Landau Lagrangian. In order 
to describe the macroscopic current-charge fluctuations 
we must also include a time dependence in the super- 
conductor order parameter. This is not an obvious task 
but it can be done in two special limits: close to T = 0, 
and near the critical temperature 27 . In the later, the 
validity of the time-dependent Gorkov equations expan- 
sion in the gap, A s , require frequencies higher than the 
binding energy of the Cooper pairs, fku A s . In this 
limit the fluctuations have enough energy to break the 
Cooper pairs and convert them into single particle exci- 
tations, leading to a diffusive regime. As we mentioned 
above, we will work in the opposite limit T — where 
the hydrodynamic description at hu> < A s is rigorously 
valid. We introduce the low-temperature time-dependent 
Ginzburg-Landau Lagrangia n 27 ; 28 , 
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which is specially suitable for type II superconductors, 
where the penetration length is much larger than the co- 
herence length (in IjlUI) (j> is the electrostatic potential 
and all the other symbols are standard^). In the clean 
limit, v a is typically of the order of the Fermi velocity vf- 
Under the assumption that the fields vary much slower 
than the coherence length, we assume that the magni- 
tude of the superconductor order parameter is constant 
and therefore all the fluctuations come from the super- 
conducting phase, ip(t,x). Thus, the order parameter is 
written as: 



(ii) 



Despite the gauge invariance of (|1(J|) , we conveniently as- 
sume the transverse gauge, V • A = 0, in what follows. 
Notice, however, that Cn + C-gl is not complete because 
the fields and the supercurrents remain decoupled. This 
problem is solved by assuming the validity of the non- 
homogeneous Maxwell equations. The Poisson equation 
is clearly a constraint between the density of supercon- 
ducting electrons and the electrostatic potential, 



47rp s = , 



and has to be enforced using a Lagrange multiplier A(x fJ ') 
in final Lagrangian. 

In the semi-classical regime, h — > 0, the behavior of the 
fields comes from the minimal action principle 

SS = sj d 4 x {C N +C GL +Ag) = sJ d 4 xC = 0, 

with respect to the variables A^, p s , X % , ip and A. By 
minimizing with respect to A M and A, we obtain the two 
non-homogeneous Maxwell's equations, 
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c 

with the supercurrent j£ given by 

*>T,2 



(13) 



p s = -^r(hcdo<p + e*0)- 6V -X-V 2 A (14) 
m*vi 



HVep- —A 

c 



cSdoX. 
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For simplicity we will assume that A 1 • = S5ij is a diago- 
nal and isotropic tensor. The minimization with respect 
to the other fields completes the set of equations: 



1 



— n {Lj- 2 d? Ps + Ps ) = a 



vV-^(0°)V-fi3V = o 

vi hvi 



k (<9 2 X - v 2 ph V 2 X) - (5(V0 + d°A) = . 



(16) 



(17) 



(18) 



Substituting the supercurrent in the continuity equation, 
J 1 — 0, we find that 



c d p s + d l j a 



c9 V 2 A = 0. 



471-ps = 



(12) 



and therefore, in the absence of external currents A is 
either a function of space or of time, but not of both. If 
A depends on time, then p is a stationary time depen- 
dent homogeneous (or non-periodic) distribution because 
of H16| l , what obviously violates the charge conservation. 
Therefore, we conclude that 8qK = 0. 

The Lagrange multiplier defines two distinct classes of 
solutions for the charge density: (i) screening-like when 
A(x) ^ 0, and (ii) plasmon-like when A = 0. In the 
first case, p(u,k) oc /(k)J(tj), where /(k) depends on 
the boundary conditions, has no dynamics and describes 
the physical response to a boundary perturbation such 
as a static squeeze of the crystal or the introduction of 
a charge probe. In case (ii), the electrons can oscillate 
freely with the plasmon frequency, allowing the existence 
of normal modes. 



IV. COLLECTIVE MODES 

We consider the problem of a layered solid (such as 
the TMD) made by an infinite number of weakly inter- 
acting planes. In the continuum limit this problem be- 
comes an effective 3D model with spatial anisotropy in 
the direction perpendicular to the planes, say the z axis. 
It is convenient to define the diagonal anisotropy ten- 
sor r* = (1, 1,t), where r is the anisotropy parameter, 
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(we introduce the notation d % = and A 1 = T l jA>) 

and rewrite the Ginzburg-Landau Lagrangian as a time- 
dependent Lawrence-Doniach modeli^: 



-a|*| 2 - 
1 

2m*v 2 
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(19) 



The model leads to the anisotropic version of equation 



r(9°)V 



hv 2 



5°</» - 



■V-A = 0, (20) 



and to the appearance of Josephson currents between the 
planes, represented by an anisotropic supercurrent: 



m* 



e ~ 

ft V<^,x) A(i,x) 

c 



c<59 X(t,x). 

(21) 

From now on, we use an overhead symbol to represent 
the in-plane vectors x, q, X,V and etc. In our notation, 
we define q as the in-plane component of the momentum 
k, such that k = (q, k z ) and k = (q, rk z ) . 

In highly anisotropic layered compounds the absence 
of piezoelectricity along the z axis is justified by the weak 
distortions of the ions in this direction. The z-component 
of the sound velocity is also very small due to the weak 
elastic coupling between planes. In the v p h z /v v h — * 
limit, the phonon dispersion is uj p h = v p hq and the in- 
plane phonon equation (|18f) has a cylindrical symmetry 



J ph) 

iS 



q<p(w,k) + -A(u,k) 



= 0. 



(22) 



Given the plasmon frequency, u> p (k), the equation for 
the charge density in Fourier space becomes: 

^[-uj 2 +w{j(k)] p.(u,k) = 8 7 r 2 o;2( k )A(k) ( 5(w). (23) 

The screening-like solution is: 



p s (u,k)=8n*^ 



5(u) 



(24) 



for a non-zero A(k). After a straightforward calcula- 
tion, the simultaneous solution of equations (|12H — (I14II . 
P ty -ffi ty and 121 yields p s {ui,k) = A (k) S {k 2 + fc 2 ), 
for some Aq function, with the characteristic momentum 



the piezoelectric coupling is larger than a critical value 
(k 2 < 0), the screening is suppressed and a quasi-static 
charge modulation should be observed, in such way to 
minimize the clastic energy. This analysis is confirmed 
by introducing an external charge probe Q at the origin: 
we re-obtain the well known Thomas- Fermi result*^, 



Pa(v, k) 



2tt 



(26) 



which gives a screened potential for A; 2 > and recovers 
the metallic case when the cut-off of the in-plane band 
dispersion s (of the order of inverse of the in-plane lattice 
spacing) is taken to infinity. 

For a finite cut-off s we have various cases. When 
fco > s, the screening is limited to the direction perpen- 
dicular to the planes. When k 2 < and |fc | < s, the 
system does not show any screening since the potential 
decays like cos(\k \r)/r for large r. The last possible 
case, feg < and |fco| > s is not physical in this theory 
and gives a purely imaginary response. We verify that 
the phase of the superconducting order parameter is free 
of fluctuations and satisfies the zero-field vortex equation 
(y) 2 ip = in the whole A^O class. 

The collective modes follow from a slightly different 
calculation. We start from the plasmon solution of (|23|) , 



p s (u),k) = p (k)5(u - w p (k)) . 



(27) 



where /?o(k) is a function which depends not only on the 
boundary conditions but also on the initial conditions 
driven by the perturbation. In the non-relativistic limit 
vf/c — * 0, the field solutions in the momentum space are 
(the details of the calculation are given in the Appendix) : 



p s (u, k) = p(q) 5(q - q ) 5{k z ) 5{u - w p (k)) . 



(u, k) = An 



Ps{u,k) 
k 2 



tp(u>, k) = 4iri 



ue* p s (a>,k) 1 

H k 2 W 2 _ v 2 ~ k 2 



5 q p s (u,k) 



(28) 



(29) 



(30) 



(31) 



kl = 4. {^3 - J± V (25) A( W| k) = 4. y-^l ^ k) , (32) 



KV P h 



where we have labeled 



The fcg > case leads to exponentially decaying solutions 

that describe the screening induced by proper boundary X>(u> k) = 1 — It * 3 * ^° \ l 4^ 



conditions like the squeeze or shear of the crystal. If m* a; 2 — v 2 k 2 K to' 2 — v 2 h q 2 
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Equation (|28|l represents the bulk plasmon mode k z = 0, 
which prevails over the rest of the plasmon band. This 
is in agreement with a general result valid for a stack 
of layers coupled by Coulomb interactions and by coher- 
ent hopping terms between adjacent planes. The pres- 
ence of inter-layer charge transfer induces a dimensional 
crossover to a 3D system, which dominates the long wave- 
length spectrum. As we show in the end of the Appendix, 
this effect is conditioned to the assumption that t^O. 
The inverse of the in-plane modulation scale q$ is given 
by the equation T>(lu p , q ) — 0. 

Noting that at T = the electronic density n is twice 
the density of Cooper pairs 'J'o, we realize that the quan- 
tity 47re* 2 5'o/ m * can be conveniently written in the form 
of the expression that gives the square of the plasma fre- 
quency, CI 2 = An e 2 n/m. In general grounds, n is given 
by the sum rulei^ 



2 f°° e 2 n 
— / duj to Im e(ui, k) — An /(<z/<z) 

7T J TO 



(33) 

where Im e is the imaginary part of the electronic sus- 
ceptibility, which is given in RPA by Ime(cj,k) = 
— Vo(q, fc z )Imn°((x>, q, ju,) and f(q/q) is the anisotropy 
function due to the shape of the Fermi surface, with 
f(q/q) = 1 in the isotropic case. For an anisotropic nodal 
liquid with a small pocket we find that f(q/q) = <j 2 /q 2 , 
as it may be easily checked by replacing the leading intra- 
band polarization function^ 



into (J2SJ, giving 



g/k* F <t:i kX 



nhv A d y/v 2 F q 2 



Air 



e* 2 ^ 2 

TO* 



u , 



with ujq = ^2kpVpe 2 / (de' hvA) , where e' is the back- 
ground susceptibility renormalized by the small inter- 
band termsii To simplify the analysis we consider the 
isotropic case only, where the rigid approximation is 
more rigorous. Using the k z = bulk plasmon mode 
uj 2 = uj 2 + (voq) 2 , we find that 
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where a is a sign function, b = (g - 



A ^ 

9 = An — 

K 



(34) 



])} and 



(35) 



is the effective piezoelectric coupling. 

Noting that v s > v p h , because the ions are much slower 
than the electrons, we identify two distinct cases, (1) v > 
v s and (2) vq < v s . The normal modes are clearly absent 
in the regime (1) for any value of g. On the other hand, 



these modes are possible in the regime (2) for non zero g. 
Naturally, they would not be physical for g < (v 2 h /v 2 )uj 2 , 
where the interactions are screened. 

The clean limit of the superconductor gives exactly 
v s = vp/y/3. In the particular case of an isotropic nodal 
liquid with a small pocket, the intra-band bulk plasmon 
disperses with vq = (v / 3/2)i>i^ and therefore it cor- 
responds to the regime (1) mentioned above. This is 
also the case of the metal (vq — \/3/5vf ), and of doped 
graphite (vo — (V3/2)vp). Therefore, in the particular 
case of TMD, we should not observe any essential differ- 
ences in the plasmon modes in comparison to the LEG 
(where g = 0) because of the piezoelectricity. 



A. Experimental results 

There is a large amount of experimental literature deal- 
ing with the observation of commensurate (CDW) or in- 
commensurate charge order in TMD. Neutron and x-ray 
diffraction measurements in TaSe2 and NbSe2 reveal the 
existence of Bragg peaks at incommensurate wave vec- 
tors Qi = (1 — 5i)bi/3, where hi (|b»| = An/{\/ia)) 
are the three reciprocal vectors with hexagonal sym- 
metry, a is the lattice spacing and Si < 0.02 is the 
incommensurability 5,9 . This state is called a triple CDW 
phase. 

In TaSe2, the phase diagram temperature vs. pres- 
sure, P, is very rich^S with three different phases: (1) 
a high temperature hexagonal incommensurate CDW 
phase (HCDW) where the three ordering vectors have 
Si ^ (i = 1,2,3); (2) an incommensurate stripe phase 
where one ordering vector is incommensurate, say Si =/= 0, 
but the other two are commensurate, (52,3 — 0; (3) a 
commensurate CDW phase (CCDW) where Si — (i — 
1,2,3). The transition from the undistorted phase (nor- 
mal) to HCDW occurs at a pressure independent temper- 
ature, Tn_h ~ 120JC, up to pressures of 4.5 GPa; the 
transition between HCDW to stripe phase Th-s ~ 110K 
is also roughly pressure independent; the transition be- 
tween stripe phase and CCDW, Ts-c(P), is highly pres- 
sure dependent and vanishes at P — P c ps 1.8 GPa. Thus, 
there is a quantum phase transition (T = 0), as a func- 
tion of pressure at P = P c . By further application of 
pressure there is a reentrant CCDW phase that has not 
been fully studied and will be not discussed here. We will 
concentrate on the nature of the quantum critical point 
(QCP) at P = P c . 

The nature of the stripe phase has been discussed 
by McMillan^ , and others^ from the phenomenologi- 
cal point of view as a result of the formation of topo- 
logical defects of the complex CDW order parameter. 
This stripe phase can be easily observed with electron 
microscopj*^. The microscopic nature of this phase is 
still unknown but our results here suggest that it may 
have its origins on the piezoelectric coupling in the solid. 
Indeed, we have found that there is a charge modu- 
lation in these materials when k 2 < (see eq.(26)). 
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This modulation can be thought as an incommensura- 
tion Si = 3 y/— fcp / 1 b j | . Equivalcntly, using definition H35|) 
there is a critical coupling constant g c so that an extra 
charge modulation appears (Si > 0) when k$ < or 



g > g c = 4?r 



i) 2 c*2,t,2 



(36) 



Recall that TaSe 2 (T c » 0.1*0 and NbSe 2 (T c w 8.3*0 
are both superconductors^ whose T c increase under ap- 
plication of pressure^ and therefore f ^ at 7 1 = 0. 
Notice that both the piezoelectric coupling, S, as well as 
the sound velocity, v p h , are also monotonically increasing 
functions of pressure. If we assume that in TaSe 2 we have 
g < g c at T = and ambient pressure, then the quantum 
phase transition can occur as a function of pressure as 
long as g/Vp h is an increasing function of pressure. In 
NbSe 2 , however, the system is always incommensurate 
indicating that g > g c even at ambient pressure. Since 
the T = phases seem to be directly connected with 
the stripe phases that are observed at finite temperature 
we can immediately conclude that the existence of these 
stripe phases have to do with the piezoelectric coupling 
in these materials. Besides, recent experiments report 
a significant electrostatic modulation of T c in epitaxial 
bilayers composed by a HT C cuprate and a polarized in- 
sulator deposited on it24. The connections between TMD 
and HT C are remarkable and we believe that this experi- 
mental result could be numerically calculated by relaxing 
the rigidity in the amplitude \& of the superconductor 
order parameter. In addition, we suggest that similar 
electrostatic devices combined with neutron scattering 
could test experimentally the role of the piezoelectricity 
in TMD stripe phases. 

Neutron scattering measurements have shown the soft- 
ening of the phonon optical mode Si in NbSe 2 and TaSe 2 
at the Qi position for a wide range of temperatures?. 
This behavior was interpreted as a coupling of the opti- 
cal phonons with the charge order. Constant-Q scans at 
Qi have shown that the phonon energy gap in both cases 
is of the order of 10 meV— . The coupling of the plasmons 
with the lattice vibrations in the present theory drives 
the long wavelength phonons to lock their frequency with 
the frequency of the k z — bulk plasmon. In compar- 
ison to the metallic case where the plasmon modes are 
rigid (in the sense that their wavelength has exactly the 
size of the system), we have found that the increase of 
the elastic energy due to the piezoelectric coupling may 
give rise to elastic plasmons which oscillate in resonance 
with the optical phonon mode. Despite not observable 
in TMD, we believe that this effect would be the macro- 
scopic manifestation of the plasmon-phonon resonance 
observed experimentally in these materials. 



CONCLUSIONS 



calculation in the continuum limit of a highly anisotropic 
material. This procedure is analogous to the Lawrence- 
Doniach effective model for an infinite stack of layers. 
Despite the evident interest in TMD, the calculation is 
sufficiently general and be could be applied to any super- 
conducting layered compound at zero temperature with 
broken lattice inversion symmetry. 

We have demonstrated from the elcctrodynamic point 
of view that superconductivity and piezoelectricity can 
coexist. Metallic screening is observed when the effective 
piezoelectric coupling g is smaller than a critical value 
g c , with the Thomas-Fermi momentum reduced by the 
increase of g. Above the critical coupling, the system 
is not screened and a long-range charge modulation is 
expected to appear as the response to a local quasi-static 
charge unbalance, which can be created by squeezing the 
crystal. We have shown that piezoelectricity is possibly 
related to the mechanism behind the stripe formation in 
TMD. Besides, we have investigated the existence of zero 
temperature normal modes arising in the presence of low 
energy bulk plasmons, which dominate the spectrum of 
collective excitations. We have also conjectured that the 
k z = plasmon mode behind these excitations comes 
from the contribution of the intra-band excitations of a 
pocket opened around the nodes of TMD by low energy 
coherent hopping terms between adjacent planes. These 
pockets could be also generated by doping TMD with 
intercalating materials. 



Acknowledgments 

A. H. C. N. thanks C. Varma for stimulating this 
work. We thank A. J. Leggett for many discussions on 
the problem of screening in metals and superconductors. 
B.U. and G.G.C. thank E. Miranda and Y. Copelevich 
for many helpfull discussions. They also acknowledge 
Fundagao de Amparo a Pesquisa do Estado de Sao Paulo 
(FAPESP), Brazil, project number 00/06881-9, for finan- 
cial support. B.U. thanks the Dept. of Physics at Boston 
University for the hospitality. 



APPENDIX 

Here we derive in detail the calculation of equations 
— (|32[) . Starting from the plasmon- like solution 



Ps{u>, k) = p Q (k) S(u - w p (k)) . 



(A.l) 



for a general function po(k) and replacing it in the the 
Poisson equation (JT5J, we get 



We have investigated the collective excitations of the 
electrons in a layered superconductor by a semi-classical 



(u;,k)=47r^^£(u;-w p (k)). (A.2) 
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Next we separate the ^-component of (|13|1 and 1211) . 

-fc a ^(w,k)+ — -fc 2 -— t A (w,k) 

c \ cr c 2 to* y 

H -i r A: z fi k) = , 

c m* 



and the <p equation H2UI) 



-k 



<p(co, k) 



(w,k)+k- A(w,k) =0 



Noting that k • A = k z (r 2 - 1 ) A z , for non-identically zero 
po(k) we find 

At n An Lok z p (k) 



1-47T- 



c /c 2 .F(u;,k) 



to*i> 2 Q(uj,k) 



c 2 jF(w,k) 



(A.3) 



and 



l-(r 2 -l) 



c 2 ^(w, k) 



(A.4) 



where we have defined: 



, 2 47T e* 2 ^§ 



TO 



and 



- p* 2 * 2 i- 2 



to*c 2 !F(u>, k) 



r(r 2 -l) 



The remaining results are derived from the phonon 
equation QZ} 

k) = -i - 1 2 (q4>(u, k) + -A(u;, k)) 

and from the combination of the in-plane components of 
(PI) and J2H), 



4ti- e* 2 ^l 



_ 

C 2 C 2 TO* 



A(w,k) + -g0(w,k) 
c 



47T e*fn 47T 

H -i qhcp(uj, k) H i w J k) = . 

C TO* C 

After a straightforward calculation, we encounter that 



X(0,,k) = -47TE--I ^ i( U -^(k)) 

K A: 2 (lo 2 - uj z ph ) 
u 2 D(w,k) 



1 - 



c 2 £(w,k) 



(A.6) 



where we have labeled 



£>(w, k) = 1 + 4vr 



S 2 



47T 



K LO — LO 



ph 



l-(r 2 -l) 



e* 2 ^ 2 1 
m*v 2 Q(u>,k) 

2 u2 

c ™ 5 



c 2 .F(w,k) 



£(w,k) = — -fc 2 - 



47re* 2 * 2 , 4tt ((5cj) 5 



c 2 k (w 2 - UJ 2 h ) 



Next, we substitute these results in the superconduct- 
ing charge density definition (|14|) 



Po(w,k) 



[Kilo Lp(u>, k) + e* 4>(lo 1 k)] 



-M<5<f-X(w,k) =0, 



what yields 



p (k)S(u;~Lu p (k))\l- 

I,-' 2 



4tt e* 2 * 2 , 



v 2 g(u,k) 



x l-(r 2 -l) 



c 2 ^(w.k) 



4tt- 



5 2 1 



1 



- 1 



u) 2 V{u,k) 



k k 2 (lo 2 — Lo 2 h ) l c 2 £(lo, k) 
= (A.7) 

Applying the transverse gauge to (|A.3(I and (|A.5(I . 

47T LO 



k- A 



k 2 



po(k) 5(u> - Lo p (k)) 



k 2 



F{w, k) 



1 - 4tp 



e* 2 * 2 



m*v 2 Q(u>,k) 



x 1 - (r 2 - 1) 



c 2 k) 



= 



2 Z>(«;,k) 

(A.8) 



This way, we conclude that /Oo(k) is on the form 

Po(k) = p{q)S(q - q )S(k z - k z0 ) , 

where qo and fc z q are the zeroes of l|A.8ll and (|A.7|I . In the 
non-relativistic limit vf/c — > 0, !F(u> p ,k) ~ £(u> p ,k) — > 
— fc 2 and 

^(^,k)->^-fc 2 



r, n 4?r ^po(k) D(w,k) 



(A.5) 



p* 2 * 2 i 

2?(w,k) -> 1-47T— 



■ 47T- 



1 



TO* CJ 2 -V 2 k 2 K LO Z - LO Z ph 
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In this limit, 1A.8|) and (|A. 7|) simplify respectively to Comparing both and recalling that k 2 = q 2 + r 2 k 2 , then 

r (1 — t) k 2 — 0. For t ^ 0,1, we immediately see that 



S 2 a 2 e* 2 *n Q 2 + rk 2 \ k z0 — 0. Substituting this result in one of the expressions 



'"' " j2_ U! 2^ ^ m * ^2 ^2^.2 y above and integrating in ui, we find that V(aj p , go) = 0. 

xp s (o;,k)<5(u;-u; p (k)) =0 (A.9) 



5 2 q 2 e* 2 V 2 k 2 \ 
+ 4tt ^ — 5 4tt =7- 

k a;^ - m* w 2 _ ^2 y 

X/5 s (ct>,k)<S(w-u; p (k)) = 0. (A.10) 
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